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Abstract 



. The absorption properties of the first low-mass protogalaxies (mini-halos) forming at high 



redshifts in the 21-cm line of atomic hydrogen are considered. The absorption properties of 
these protogalaxies are shown to depend strongly on both their mass and evolutionary status. 
The optical depths in the line reach ^0.1-0.2 for small impact parameters of the line of sight. 
When a protogalaxy being compressed, the influence of gas accretion can be seen manifested in a 



non-monotonic frequency dependence of the optical depth. The absorption characteristics in the 
+3 . 21-cm line are determined by the thermal and dynamical evolution of the gas in protogalaxies. 

Since the theoretical line width in the observer's reference frame is 1-6 kHz and the expected 
separation between lines 8.4 kHz, the lines from low mass protogalaxies can be resolved using 
ongoing and future low frequency interferometers. 

>: 

q ; 1 Introduction 

One of the means of studying the history of intergalactic space before the epoch of secondary 
f*"*- 1 ionization (reionization) of hydrogen is radio observations in the 21-cm line of neutral hydrogen [1- 

3]. At least three methods for studying this epoch have been proposed: (1) observations in emission 
or absorption in the neutral intergalactic medium against the global cosmic microwave background 
(CMB) signal [4, 5], (2) statistical studies of the angular distribution of the signal fluctuations [6, 
7], and (3) observations of absorption in the spectrum of a distant strong radio source, which may 
reveal a "forest" in the 21- cm line (analogous to the La forest in the optical spectra of quasars 
and galaxies) [8-10]. The first two methods involve mapping spatial (position on the sky) and 
temporal (redshift) properties of the intergalactic medium, while the latter provides information 
only along a line of sight as a function of redshift. Since the spatial scales corresponding to the 
angular resolutions of ongoing or futture radio telescopes, such as LOFAR0 are roughly 1 Mpc, only 
large-scale structures, such as giant ionization regions around massive galaxies and quasars, can 
be observed. Measurements of absorption at different redshifts can be used to resolve small-scale 
structures, individual protogalaxies (mini-halos), and ionization zones in low-mass protogalaxies. 

The effect of external ionizing radiation on the properties of the 21-cm forest was studied in 
[9-14]. Mainly, the dark-matter profiles in protogalaxies have been taken to be steep (cusped) 
[15]. However, observations of dwarf galaxies in the local Universe indicate a "flat" distribution of 
dark matter (a truncated isothermal sphere with a flat core) [16]. Recent theoretical studies have 
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also found that the first protogalaxies could have flat dark-matter profiles [17-19]. The first halos 
were assumed in [9, 10] to be virialized and isothermal, and their absorption characteristics were 
calculated without accounting of the gas and dark matter evolution. Some evolutionary effects 
were included in [14], such as a possible influence of gas accretion and star formation, although the 
gas and dark matter profiles of minihalos before the first stars have appeared were assumed to be 
time-independent. It is well known though that the density inside the perturbation from which a 
protogalaxy does form increases by more than a factor of 200 in the non-linear stage, facilitating 
cooling of the gas and a decrease of its ionization fraction [20-23] . It is obvious that such appreciable 
evolutionary variations in the physical properties of the gas could substantially affect absorption 
properties of the first protogalaxies. With this in mind, Meiksin [13] considered possible variations 
in the absorption properties caused by the virialization process of the first protogalaxies. 

Here, we investigate the absorption properties of the first protogalaxies within a self-consistent, 
one-dimensional, spherically symmetric description of their thermal and chemical evolution. In the 
heirarchical scenario of galaxy formation, the first protogalaxies able collapse and lead to formation 
of the first stars had low masses. The minimum mass of the first protogalaxies, ~ 10 7 Mq 
at redshifts z ~ 10, was apparently determined by molecular-hydrogen cooling [23]. We will 
study the 21-cm profiles for mini-halos with masses near this minimum mass, with the aim of 
using possible differences in the 21-cm absorption characteristics as a tool to answer the question 
of whether the minimum masses of halos in which star formation began were indeed close to 
M h ~ 10 7 M Q at z ~ 10. Section 2 describes the model used, Sections 3 and 4 present our results 
and conclusions. The computations assume a ACDM cosmological model with (Qq, f2 m , £lb,h) = 
(1.0, 0.76, 0.24, 0.041, 0.73), and a relative number density of deuterium of n[D)/n = 2.78 • 10~ 5 



2 Evolution of minihalos 

We modeled evolution of the first protogalaxies using a one-dimensional Lagrangian approach 
similar to that proposed in [23]. It is difficult to fully take into account mergers of protogalaxies 
and other essentially non-one-dimensional dynamical effects in such a simple model. This would 
require self-consistent, three-dimensional cosmological simulations with high resolution, which is 
not currently possible. By necessity therefore we will use simplified model approximations that are 
commonly applied to reduce multi-dimensional problems to spherically symmetric ones. 

2.1 Dynamics of Dark Matter 

The description of dark matter in our model is analogous to that of [23]. It is assumed that 
dark matter with the mass Mdm = ^DAiMhaio/^m is contained inside the radius Rt r and has a 
truncated isothermal sphere profile with a flat core of radius R CO re, here £l m is the matter density 
and Qdm = &m — ^b- The ratio £ = R core /R v i r is assumed through the paper to be 0.1; R V i r 
is the virial radius of the halo, outside the radius Rt r , the dark- matter density is equal to the 
background cosmological value, pdm = Po^da/(1 + z) 3 , where po = 1.88 x 10~ 29 h 2 is the critical 
density; the radii Rt r and R cor e evolve in time [23]. Such a description is usually accepted to 
imitate the evolution of density perturbations (see, e.g., [25]). Before the "turnaround" point, i.e. 
the point where the perturbation separates from the Hubble flow, the density within the truncated 
sphere, R < Rt r , varies with time similarly to what approximated in [21] 



where A = (1 + z v i r )/{\ + z). We accept in ([I]) the numerical coefficients in the exponent slightly 
different from the ones used in [21], in order to achive a better agreement between (fTJ) and the 
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exact solution in the vicinity of the turnaround point. The approximation (1) violates near the 
virialization point, however this is not important for our purposes, since it is not used beyond the 
turnaround point. 



2.2 Gas Dynamics and Chemical Kinetics 

For modeling evolution of the gas (i.e., the baryonic component of a protogalaxy) a one-dimensional 
spherically symmetric Lagrangian scheme similar to that described in [26] is used. The standard 
resolution in the computations is 1000 cells (spherical shells), which provides a sufficiently good 
convergence of the scheme. 

Chemical kinetics of the primordial gas includes the following main components H, H + , e, H - , 
He, He + , He ++ , H2, H^, D, D + , D _ , HD, HD + . The corresponding reaction rates are taken from 
[27, 28]. In the energy equation radiative losses typical for the primordial gas: cooling with recom- 
bination, collisional excitation of hydrogen and helium, free-free transitions, Compton interactions 
with CMB photons [29], and molecular cooling via H2 [28] and HD [30, 31] arc accounted. 

All models begin at the redshift z = 100, when the density perturbations are small and obviously 
linear. The initial parameters for z = 100, i.e., the gas temperature, ion and molecular composition, 
etc, are found from a one-zone calculateion starting at z = 1000 with the standard gas parameters 
for the end of the recombination epoch: T gas = TqmBi x\H] = 0.9328, x[H + ] = 0.0672, x[D] = 

2.3 x 10~ 5 ,x[D + ] = 1.68 x 10~ 6 (see [23, Table 2] for details and references). The initial size of the 
computational domain is 10 virial radii for a protogalaxy of a given mass. 

2.3 The 21 cm Line 

Only collisional excitation by atoms and electrons is taken into account when computing the spin 
temperature for the 21-cm line [32]: 

T _ Tqmb + VcTk , x 

s ~ i + y c {) 

where Tcmb is the temperature of the CMB, y c is a function determined by the efficiency of 
collisional excitation 

Vc = ~A — T^' W 
/i -10 1 k 

T* = 0.0682 K is the energy for the transition between sublevels of the hyperfine structure, 
A10 = 2.87 x 10~ 15 s" 1 is the transition probability, C10 = kionH + 7e^e is the rate of colli- 
sional de-excitation by hydrogen atoms and electrons (we neglect the contribution from protons), 
the approximations for k\o of [33] and for j e from [34] are used. The computations do not take 
into account the influence of resonance ultraviolet radiation (UV pumping) on the level populations 
[35], since at redshifts z = 10 — 15 this effect can be important only locally near galaxies with star 
formation. At lower redshidt the regions of penetration of UV radiation from the stellar population 
become substantially greater and overlap and manifest reionization of the Universe. Moreover, in 
order for the contribution from UV pumping to the variations of the spin temperature to be com- 
parable to the contribution from collisions, the flux of resonance UV radiation must be fairly high: 
Fhya ~ 0.1 erg s _1 cm~ 2 (this estimate is for typical temperatures T ~ 200 K and gas densities 
n ~ 1 cm~ 3 in protogalaxies at z ~ 10). 

3 Absorption in the 21 cm line 

We will study two observable quantities that describe the absorption properties of the gas in 
protogalaxies: the optical depth at frequency v along the line of sight, and the total absorption 
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Figure 1: Number density, temperature, and velocity (from top to bottom) of the gas in proto- 
galaxies with masses of 10 5 , 10 6 and 10 7 M Q (from left to right) at times zt a = 15.5 (dot-dashed 
curve), z = 12 (dashed curve), and z v = 10 (solid curve). 

within the protogalaxy, i.e., the line equivalent width. We consider the properties of low-mass 
protogalaxies (mini-halos) with masses of 10 5 , 10 6 and 1O 7 M0, which separated from the Hubble 
flow at redshift zt a — 15.5 and virialized by z v c± 10. 

Evolution of the first protogalaxies depends critically on their masses: low-mass protogalaxies 
cannot cool rapidly and form stars after their virialization. With increasing protogalaxy mass the 
number density of hydrogen molecules in the central regions grows, facilitating cooling of the gas 
and formation of the first stars (see [20-23, 36] for more detail). Figure 1 presents radial profiles 
of the gas density, temperature, and velocity in protogalaxies with masses of 10 5 , 10 6 and 1O 7 M 
at times close to the turnaround [zt a = 15.5) and the virialization [z v = 10), and also at a time 
corresponding to the intermediate redshift [z = 12). When virialized at z v = 10, only protogalaxies 
with masses M > 1O 7 M can form stars [22, 23]. In the process of virialization, gas density in 
the central regions already reaches n £ 10 8 _3 at z ~ 11, which essentially corresponds to the 
collapse and the onset of formation of the first pre-stellar cores. Therefore, for the virial state of 
the mass M = 1O 7 M we present the absorption characteristics corresponding to z ~ 11, while for 
the two other masses the characteristics are calculated for z v = 10. As noted above, in lower-mass 
protogalaxies, the prestellar cores either form after virialization or do not form at all. 

The optical depth at frequency v along a line of sight passing a distance a from the center of 
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Figure 2: Optical depth for protogalaxies with masses of 10 5 , 10 6 and 10 7 M Q (from left to right) 
at times zt a = 15.5 (upper panel), z = 12 (middle panel), and z v = 10 (lower panel). The curves 
correspond to the line of sight impact parameters a = 0.1, 0.3, 1, 1.5, 3r v i r (from right to left), and 
the thin dashed line to the background optical depth. 

the protogalaxy is 

_ 3h p c 3 A 10 [°° ]jd n H i{r) 
Tu ~ 32ttH 2 J-oo d V& 2 (r)T s (r) 6XP 

where r 2 = a 2 + R 2 , v(v) = c(y — vo)/i>o, vi(a, R) is the gas infall velocity projected on to the line 
of sight, b 2 = 2kTk(r) / m p is the Doppler parameter, and the remaining notations are standard. 

Figure 2 shows the optical depths for lines of sight passing at several impact parameters for 
protogalaxies with masses of 10 5 , 10 6 and 10 7 M Q at times close to the turnaround (zt a = 15.5) 
and the virialization time (z v = 10), as well as at a time corresponding to the intermediate redshift 
z = 12. Considerable differences in the line absorptions are clearly visible, both for protogalaxies 
with different masses and during the evolution of a protogalaxy of a given mass. At z = 15.5 (upper 
panel in Fig. 2), the gas densities and temperatures in the protogalaxies with different masses are 
nearly the same (Fig. 1), and the differences in the line profiles are determined by the distributions 
of the gas velocity: with increasing protogalaxy mass the accretion velocity increases from ~0.2 
km/s for 1O 5 M to ~2 km/s for 1O 7 M . The small peak in the optical depth near u — uq ~ 40 kHz 
and the drop to zero at higher frequencies for M = 1O 5 M is, from one side, due to the drop in 
the temperature at the periphery of the protogalaxy and the transition to the background value 
[a decrease in the denominator in the exponent (4)], and from the other due to the transition of 



[v(y) -v l {a,R)f 
b 2 {r) 



(4) 
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the velocity profile to the Hubble expansion [an increase in the numerator in the exponent (4)]. In 
more massive protogalaxies, a similar peak in optical depth is located at higher frequencies. 

The absorption characteristics change considerably as the protogalaxy is compressed. At red- 
shift z = 12 (middle panel in Fig. 2), the absorption lines become broader, and features due 
to the gas dynamics during the collapse of the protogalaxies appear. In low-mass protogalaxies 
(M = 10 5 Mq), the enhancement of gas velocity in the accretion shock is modest, so that the ab- 
sorption line profiles for the various impact parameters of the line of sight broaden only slightly. In 
higher-mass protogalaxies (M ^ 1O 6 M ), the accretion wave is stronger, and a stronger jump in 
the temperature and density arises, which manifests a nonmonotonic frequency dependence of the 
optical depth for the lines of sight intersecting this region (a £ r v i r ) as, more precisely, a decrease 
in the optical depth in the line center. 

At z = 10 (at z ~ 11 for M = 1O 7 M0), the absorption-line profiles depend substantially on 
differences in the dynamical and thermal evolution of the protogalaxies (lower panels in Figs. 1 and 
2): the lines become appreciably broader, the strongest absorption is observed for a £ 0.3r\,j r , the 
influence of gas accretion is clearly seen in the protogalaxy with mass M = 10 Mq resulting in pro- 
nounced peaks in the frequency dependence of the optical depth. In more massive protogalaxies gas 
cools more efficiently, and, moreover, the shocks associated with virialization appear to be stronger 
in more massive protogalaxies. Note that the line profile for M = 1O 6 M0 at z = 10 has changed 
since z = 12: contrary to z = 12 at z = 10 the optical depth falls off essentially monotonically from 
the line center. This is due to a decrease in the accretion, and the appearance of a reverse shock 
acting to oppose the accretion (Fig. 1). In the case of higher-mass protogalaxies (M = 1O 7 M0), 
the accretion increases with decreasing redshift, while the non-monotonic dependence of the optical 
depth becomes stronger. 

Let us now investigate how thermal evolution affects the absorption- line profiles. For this 
purpose, we calculate the evolution of a protogalaxy with mass M = 10 6 Mq without radiation gas 
cooling in atomic and molecular processes (i.e. we consuder the adiabatic case), and compare it 
with the results described above (Fig. 3). In the presence of radiation cooling the central regions 
of the protogalaxy are seen to cool efficiently during virialization, and the radial dependences differ 
considerably from the adiabatic case, but these differences are substantial only in the inner region 
of the protogalaxy, a £ 0.5r„j r , where the effects of cooling are stronger due to the higher density. 
As a result, the absorption line profiles will also differ in this region as seen in Fig. 4. 

Another observational characteristic is the line equivalent width, (Az/) & s = (Ai>)i/(l+z), where 
the rest-frame equivalent width is 



where tigm is the background optical depth of the neutral intergalactic medium. Figure 5 presents 
the observed line equivalent width for absorption in protogalaxies with masses of 10 5 , 10 6 and 
1O 7 M0 at redshifts z = 15.5, 12, and 10. The equivalent width depends only weakly on the impact 
parameter of the line of sight at z = 15.5, but in the process of formation of a protogalaxy it grows 
at a ^ r v { r , reaching 0.1-0.4 kHz for M = 10 6 — 10 7 Mq. The line width of the perturbations after 
separation from the Hubble flow in the observer's rest frame (i.e. at z = 0) equals to 1-6 kHz. Since 
the distance between the lines is ~8 kHz at z = 10 [14], and the spectral resolution of 0.76 kHz 
can be reached with existing and planned low-frequency interferometers, such as LOFAR and the 
ska| the lines of low-mass protogalaxies might be quite resolvable with these instruments. 

Note that when describing the initial growth of a density perturbation after recombination, 
it is assumed that the velocities of the dark matter and baryons are similar and their relative 
velocities low [37]. This picture has been recently supplemented with a possibility of supersonic 
relative motions of baryons and dark matter with velocities v s {z) = v s ,i/(l + z), where v s ,% is the 

2 http://www.skatelescope.org/ 
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Figure 3: Density, temperature, and velocity (from top to bottom) in a protogalaxy with mass 
1O 6 M at times zt a = 15.5 (dot-dashed curve), z = 12 (dashed curve), and z v = 10 (solid curve), 
with (thin curves) and without (thick curves) taking into account cooling of the gas. The difference 
is clearly seen at the virialization time z v = 10; these curves coincide at the other two times 
considered. 
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Figure 4: Optical depth for a protogalaxy with mass 1O 6 M at time z v = 10. The curves 
correspond to the line of sight impact parameters a = 0.1, 0.3, 0.5, 0.8, lr V i r (from right to left) 
with (thin curves) and without (thick curves) radiation gas cooling. The difference is clearly seen 
for a curves coincide for higher values of a. 
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Figure 5: Observed line equivalent width for protogalaxies with masses of 10 5 , 10 6 and 1O 7 M 
(from top to bottom) at times zt a = 15.5 (dot-dashed curve), z = 12 (dashed curve), and z v = 10 
(solid curve). 
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rms relative velocity of the gas at the epoch of recombination [38]. It is assumed that, during the 
compression of protogalaxies, their gas moves with an effective speed v e ff = yj cj. + v'j [c s is the 
sound speed) [39] . We carried out computations of the evolution of protogalaxies with masses of 
M = 10 5 — 10 7 Mq taking into account v s (z), adopting the characteristic value v s ^ = 30 km/s 
for the epoch of recombination [38]. It was found that compared to the standard case described 
above (Fig. 2) the accounting of the relative motion of baryons and dark matter results in an 
approximately 10-30% decrease in the optical depth in the 21-cm line at the epoch of separation 
of the perturbation from the Hubble flow. For more massive protogalaxies this change is smaller 
and decreases by the epoch of virialization, while for low-mass protogalaxies remains considerable 
comprising about 10%. 

4 Conclusions 

We have used one-dimensional, spherically symmetric computations to investigate the absorp- 
tion properties of the first protogalaxies in the 21-cm line of HI at the lower mass boundary 
admitting the collapse of baryons and star formation, Mh < 10 7 M at z £ 10. We have studied 
the evolution of the 21-cm line profile and observed equivalent width during the virialization of 
these protogalaxies. 

The absorption properties of the protogalaxies depend substantially on both their mass and 
their evolutionary status. The optical depth in the line reaches ~ 0.1 — 0.2 at small line of sight 
impact parameters for protogalaxies with masses of M = 10 5 - 1O 6 M , then monotonically falls 
off to the background value. The 21-cm absorption profiles for halos with masses of M = 10 Mq 
and probably higher masses differ qualitatively from the profiles of the halos with masses below 
this lower limit: massive halos show a dip at the line center. This can be used to observationally 
determine the minimum masses of protogalaxies; i.e., of minihalos capable to cool radiatively, and 
enabling the collapse of baryons with subsequent star formation. 

The influence of gas accretion, manifested as a nonmonotonic frequency dependence of the 
optical depth, can be seen as the protogalaxies are compressed. The accretion wave is observed 
at large line-of-sight impact parameters for protogalaxies with masses M ~ 10 5 M Q , and becomes 
stronger for protogalaxies with higher masses (M > 1O 6 M ), leading to a larg e jump in the gas 
temperature and density, which gives rise to appreciable peaks in the frequency dependence of the 
optical depth for a £ r v i r . The absorption properties are mostly determined by the thermal and 
dynamical evolution of the gas in the protogalaxies. 

The equivalent width corresponding to the turnaround point in the observers rest frame (i.e., 
z = 0) depends only weakly on the line-of-sight impact parameter a, but grows with the forma- 
tion (virialization) of protogalaxies with masses of M = 10 6 — 1O 7 M , reaching 0.1-0.4 kHz for 
a ;$ 0.5r v i r . The expected distance between lines is 8.4 kHz at z ~ 10 [14], that is larger than the 
width of the 21 cm absorption lines from low mass protogalaxies in the observers rest frame, so 
that the lines can be resolved using low frequency ongoing and future interferometers. 
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